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Abstract. Wc study partition functions and thermodynamic limits for the 
Ising model on three families of finite graphs converging to infinite self-similar 
graphs. They are provided by three well-known groups realized as automor- 
phism groups of regular rooted trees: the first Grigorchuk's group of interme- 
diate growth; the iterated monodromy group of the complex polynomial — 1 
known as the "Basilica group"; and the Hanoi Towers group H^^^ closely re- 
lated to the Sierpihsky gasket. 
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1. Introduction 

1.1. The Ising modeL The famous Ising model of ferromagnetism was introduced 
by W. Lcnz in 1920, [TT], and became the subject of the PhD thesis of his student 
E. Ising. It consists of discrete variables called spins arranged on the vertices of a 
finite graph Y . Each spin can take values ±1 and only interacts with its nearest 
neighbours. Configuration of spins at two adjacent vertices i and j has energy 
Jij > if the spins have opposite values, and —Jij if the values are the same. Let 
|yert(Y')| = N, and let a = (cti, cttv) denote the configuration of spins, with 
<7i G {±1}. The total energy of the system in configuration a is then 



^(ct) = JjjCTjCTj, 



where we write i ^ j if the vertices i and j are adjacent in Y. 

The probability of a particular configuration at temperature T is given by 

P(a) = iexp(-/3i?(a)), 

where /3 is the "inverse temperature" conventionally defined as /? = l/lksT), and 
denotes the Boltzmann constant. 

As usual in statistical physics, the normalizing constant that makes the distri- 
bution above a probability measure is called the partition function: 



Z = ^exp(-/?ii;(a)). 
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One can rewrite this formula by using exp(A'(Ti(Tj) = cosh(/\) + aicrj sinli(/\), so 
as to get the so-called "high temperature expansion" : 

Z = "|~[ cosh(/3 Jj^j) + gifTj tanh{l3Jij) 

+ ^ (CT.iCTj)(criCTm) tanh(^J,j) tanh(^Ji^m) H ). 

After changing the order of summation, observe that the non-vanishing terms in Z 
are exactly those with an even number of occurrences of each ai . We can interpret 
this by saying that non-vanishing terms in this expression are in bijection with closed 
polygons of Y, i.e., subgraphs in which every vertex has even degree. Consequently 
we can rewrite Z as 

(1)Z= (ncosh(/?J,,)) -2^ n tanh(/3J,,), 

j~J X closed polygon of Y {i,j)eEdges{X) 

where in the RHS we have the generating series of closed polygons of Y with 
weighted edges, the weight of an edge {i,j) being tanh(/?Jij). 
In the case of constant J, the above expression specializes to 

Z = (cosh(/3J))l'''^^'^''^l • 2^r^'(tanh(/3J)) 

with r'^'(z) = X^riLo -^"j where A'^ is the number of closed polygons with n edges 
in Y. (In particular, the total number of closed polygons is given by r'^'(l).) 

From the physics viewpoint it is interesting to study the model when the system 
(i.e. the number of vertices in the graph) grows. One way to express this mathe- 
matically is to consider growing sequences of finite graphs converging to an infinite 
graph. If the limit 

lim , ^"g(f"\, 
n->oo \Vert{Yn)\ 

for a sequence of finite graphs K„ with partition functions Z„ exists, it is called the 
thermodynamic limit. 

In the thermodynamic limit, at some critical temperature, a phase transition 
can occur between ordered and disordered phase in the behaviour of the model. 
Existence of a phase transition depends on the graph. In his thesis in 1925, Ising 
studied the case of one-dimensional Euclidean lattice; computed the partition func- 
tions and the thermodynamic limit and showed that there is no phase transition 
[ID] . The Ising model in 'E'^ with d > 2 undergoes a phase transition. This was first 
established for d = 2 by R. Peierls. At high temperature, T > Tc, the clusters of 
vertices with equal spins grow similarly for two different types of spin, whereas for 
T < Tc the densities of the types of spin are different and the system "chooses" 
one of them. 

The infinite graphs that will be studied in this paper all have finite order TZ of 
ramification, i.e., for any connected bounded part X of the graph there exists a set 
A of at most 7?. vertices such that any infinite self-avoiding path in the graph that 
begins in X necessarily goes through A. Finite order of ramification ensures that 
the critical temperature is T=0, and there is no phase transition in the Ising model 
(see [5]). 
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Typically, an infinite lattice is viewed as the limit of an exhaustive sequence of fi- 
nite subgraphs. This is a simple example of the so-called pointed Hausdorff-Gromov 
convergence (see Proposition 11.31 below) . Another typical case of this convergence 
is that of covering graph sequences. In this paper wc will be studying the Ising 
model on families of finite graphs coming from the theory of self-similar groups 
(sec Definition 1 1 . 21 below) . and their infinite limits. Any finitely generated group of 
automorphisms of a regular rooted tree provides us with a sequence of finite graphs 
describing the action of the group on the levels of the tree. When the action is 
self-similar the sequence converges in the above sense to infinite graphs describing 
the action of the group on the boundary of the tree. The graphs that we study 
here are determined by group actions, and so their edges are labeled naturally by 
the generators of the acting group. Different weights on the edges lead to weighted 
partition functions, with J^j depending on the label of the edge (i, j). 

1.2. Groups of automorphisms of rooted regular trees. Let T be the infinite 
regular rooted tree of degree q, i.e., the rooted tree in which each vertex has q 
children. Each vertex of the n-th level of the tree can be regarded as a word of 
length n in the alphabet X = {0,1, — 1}. Moreover, one can identify the 
set X" of infinite words in X with the set dT of infinite geodesic rays starting at 
the root of T. Now let G < Aut{T) be a group acting on T by automorphisms, 
generated by a finite symmetric set of generators 5". Suppose moreover that the 
action is transitive on each level of the tree. 

Definition 1.1. The n-th Schreier graph S„ of the action of G on T, with respect 
to the generating set S, is a graph whose vertex set coincides with the set of vertices 
of the n-th level of the tree, and two vertices u, v are adjacent if and only if there 
exists s G S such that s(u) = v. If this is the case, the edge joining u and v is 
labeled by s. For any infinite ray ^ € dT, the orbital Schreier graph has vertices 
G ■ ^ and edges determined by the action of generators, as above. 

The vertices of I]„ are labeled by words of length n in X and the edges are 
labeled by elements of S. The Schreier graph is thus a regular graph of degree 
d = \S\ with g" vertices, and it is connected since the action of G is level-transitive. 

Definition 1.2. A finitely generated group G < Aut(T) is self-similar if, for all 
g G G, X Cz X , there exist h Cz G,y Cz X such that 

g{xw) ^ yh{w), 

for all finite words w in the alphabet X . 

Self-similarity implies that G can be embedded into the wreath product Sym{q) i 
G, where Sym{q) denotes the symmetric group on q elements, so that any auto- 
morphism g E G can be represented as 

g = r(go, ■ • ■ 

where r G Sym{q) describes the action of g on the first level of T and gi e G, i = 
0, (7 — 1 is the restriction of g on the full subtree of T rooted at the vertex i of 
the first level of T (observe that any such subtree is isomorphic to T). Hence, if 
X & X and w is a finite word in X, we have 



g{xw) = T{x)gx{w). 
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It is not difficult to see that the orbital Schreier graphs of a self-similar group 
are infinite and that the finite Schreier graphs form a sequence of graph 

coverings (see |12| and references therein for more information about this interesting 
class of groups, also known as automata groups.) 

Take now an infinite ray ^ e X" and denote by ^„ the n-th prefix of the word 
^. Then the sequence of rooted graphs {(!]„, ^„)} converges to the infinite rooted 
graph (E^,^) in the space of rooted graphs, in the following sense. 

Proposition 1.3. f]9], Chapter 3.) Let X he the space of connected graphs having 
a distinguished vertex called the root; X can he endowed with the following metric: 
given two rooted graphs (Yi,ui) and (12,^2)7 

Z?i.st((Yi, ui), (I2, '^2)) := inf < -;_By, (wi,r) is isomorphic to By, (v2 , I 

[r+1 - J 

where BY(v,r) is the ball of radius r in Y centered in v. Under the assumption of 
uniformly hounded degrees, X endowed with the metric Dist is a compact space. 

1.3. Plan of the paper. Our aim in this paper is to study the Ising model on the 
Schreier graphs of three key examples of self-similar groups: 

- the first Grigorchuk's group of intermediate (i.e., strictly between polynomial 
and exponential) growth (see [7] for a detailed account and further references); 

- the "Basilica" group that can be described as the iterated monodromy group 
of the complex polynomial — 1 (see [12] connections of self-similar groups to 
complex dynamics); 

- and the Hanoi Towers group iJ^^^ whose action on the ternary tree models the 
famous Hanoi Towers game on three pegs, see [8]. 

It is known [T] that the infinite Schreier graphs associated with these groups 
(and, more generally, with all groups generated by hounded automata) have finite 
order of ramification. Hence the Ising model on these graphs exhibits no phase 
transition. 

We first compute the partition functions and prove existence of thermodynamic 
limit for the model where interactions between vertices are constant: in Section [2] 
we treat the Grigorchuk's group and the Basilica group, and in Section[3]the Hanoi 
Towers group and its close relative the Sierpihski gasket are considered. 

In Section 21 we study weighted partition functions for all the graphs previously 
considered, and we find the distribution of the number of occurrences of a fixed 
weight in a random configuration. The relation between the Schreier graphs of 
H^^^ and the Sierpihski gasket is also discussed from the viewpoint of Fisher's 
theorem establishing a correspondence between the Ising model on the Sierpihski 
gasket and the dimers model on the Schreier graphs of H^^^ . 

1.4. Acknowledgments. We are grateful to R. Grigorchuk for numerous inspiring 
discussions and to S. Smirnov for useful remarks on the first version of this paper. 

2. Partition functions and thermodynamic limit for the 
Grigorchuk's group and for the Basilica group 

2.1. Grigorchuk's group. This group admits the following easy description as a 
self-similar subgroup of automorphisms of the binary tree. It is generated by the 
elements 

a — e{id,id), b — e{a,c), c = e{a,d), d=e(id,b), 
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where e and e are respectively the trivial and the non-trivial permutations in 
Sym{2). These recursive formulae allow easily to construct finite Schreier graphs 
for the action of the group on the binary tree. Here arc three first graphs in the 
sequence, with loops erased. 




c 




In general, the Schreier graph E„ has the same linear shape, with 2"~^ simple 
edges, all labeled by a, and 2""-'^ — 1 cycles of length 2. It is therefore very easy to 
compute the generating function of closed polygons of £„, for each n. 

Theorem 2.1. The generating function of closed polygons for the n-th Schreier 
graph of the Grigorchuk group is Tjl{z) = (1 + z^)^ . In particular, the number 
of all closed polygons in E„ is 2^ ~^ . 

The partition function of the Ising model is given by 

Zn ^ cosh(/3J)3-2""'"2 • 2^" ■ (1 + tanh^ (/3J))^" '"^ 
and the thermodynamic limit exists and satisfies: 

hm = I log(cosh(/? J)) + log 2 + i log (1 + tanh^ (/3 J)) . 

Proof. Is is clear that a closed polygon in S„ is the union of 2-cycles. So we can 
easily compute the number A'^^ of closed polygons with k edges in S„, for all 
= 0, 1, . . . , 2" — 2. For k odd, one has A'^^ = 0. For k even, we have to choose | 
cycles of length 2 to get a closed polygon with k edges, which implies 

So the generating function of closed polygons for I]„ is given by 

□ 

2.2. The Basilica group. The Basilica group is a self-similar group of automor- 
phisms of the binary tree generated by the elements 

a = e(6, id), b ~ e(a, id). 

The following pictures of graphs E„ for n = 1,2,3,4,5 with loops erased give an 
idea of how finite Schreier graphs of the Basilica group look like. See [4] for a 
comprehensive analysis of finite and infinite Schreier graphs of this group. Note 
also that {S„}5^]^ is an approximating sequence for the Julia set of the polynomial 
— 1, the famous "Basilica" fractal (see [12]). 
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Si •: 







In general, it follows from the recursive definition of the generators, that each E„ 
is a cactus, i.e., a union of cycles (in this example all of them are of length power of 
2) arranged in a tree-like way. The maximal length of a cycle in I]„ is 2^51 . Denote 
by aj the number of cycles of length j labeled by a in and analogously denote 
by b^j the number of cycles of length j labeled by b in E; . 

Proposition 2.2. For any n > 4 consider the Schreier graph E„ of the Basilica 
group. For each k > 1, the number of cycles of length 2^ labeled by a is 




forl<k< 
for k = Lf J 



2 



- 1 



for n odd, 
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for n even 



/or 1 < fc < f - 1 
for fc = f 

and the number of cycles of length 2^ labeled by b is 

{2"-2fc /or 1 < fc < ^ - 1 
2 fork= [f J , for n odd, 

1 fork= [f ] 

/2"-2fe /orl<fc<f-l 
o„fc = < , for n even. 

^ [2 fork=^ ■' 

Proof. The recursive formulae for the generators imply that, for each n > 3, one 
has 



^2 



= 6^-1 and = a1-^ = 2"-^ 



and in general a^^ = and 62*: — ^2 ^'^'^ particular, for each n > 4, 

the number of 2-cycles labeled by a is 2""'^ and the number of 2-cycles labeled by 
b is 2"~^. More generally, the number of cycles of length 2'° is given by 

n rfn — 2k — l in rtn — 2k 

where the last equality is true if n — 2fc + 2 > 4, i.e. for fc < ^ — 1- Finally, for n 
odd, there is only one cycle of length 2rtl labeled by b and four cycles of length 
2LtJ, two of them labeled by a and two labeled by b; for n even, there are three 
cycles of length 2 2 , two of them labeled by b and one labeled by a. □ 

Corollary 2.3. For each n > 4, the number of cycles labeled by a in the Schreier 
graph S„ of the Basilica group is 

( 2"~^+2 



and the number of b-cycles in I]„ is 




for n odd, 
for n even. 

for n odd, 
for n even. 

The total number of cycles of length > 2 is 2"^^ + 1 and the total number of edges, 
without loops, is 3 ■ 2"^^. 

The computations above lead to the following formula for the partition function 
of the Ising model on the Schreier graphs I]„ associated with the action of the 
Basilica group. 

Theorem 2.4. The partition function of the Ising model on the n-th Schreier graph 
S„ of the Basilica group is 

Zn = 2^" • cosh(^J)3-2""' . r,'=,'(tanh(/3J)), 

where T'^{z) is the generating function of closed polygons for I]„ given by 



k=l 



.2T 
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for n > 5 odd and 

fc=i 

forn>4 even. Moreover, Vf = l + z^. r^' = (1 + 22)3 andVf = (1 + 22)4(1 + ^4)^ 

Proof. Recall that Z„ = 2l^'=''*(^")l cosh(/3J)l^'*f'=''(^")l •r^'(tanh(^J)), where r^'(z) 
is the generating function of closed polygons in I]„. In our case we have |i?dges(S„)| 
= 3 • 2"-i and |yert(S]„)| = 2". 

The formulae for Tfj^iz) with n = 1, 2, 3 can be directly verified. For n > 4, we 
can use Proposition 12.21 Since the length of each cycle of I]„ is even, it is clear 
that the coefficient Al^^ is zero for every odd k. The coefficient Af^ is nonzero for 
every even k such that < fc < 3 • 2"^^. In fact. 3 • 2""^ is the total number of 
edges of (2" labeled by b and 2"~^ labeled by a). By taking the exact number 
of cycles of length 2* in E„, we get the assertion. □ 

Theorem 2.5. The thermodynamic limit lim„^oo \ v°ertf£\\ s^^'^^''- 

Proof Since |£'dges(I]„)| = 3 • 2"-^ and |yeri(S]„)| = 2", the limit reduces to 
(choosing, for example, n even) 

log(2) + I log(cosh(/?J)) + lim iHligM , 
where z = tanh(/3J) takes values between and 1. Now 

j.^^ log(r^'(z)) ^ ^.^ Et"/ 3 ■ 2"-^fe-i log(l + z2'°) + 31og(l + z^^) 

n— >-oo 2"" n— ^cxD 2^^ 



3y. l0g(l + Z^') ^ 3l0g(l + z2^) 



3 log(2) 

k=l 

giving the assertion. □ 



3. Partition functions and thermodynamic limits for the Hanoi 
Towers group H^^^ and for the Sierpinski gasket 

3.1. Hanoi Towers group ij'^'. The Hanoi Towers group is generated by 
three automorphisms of the ternary rooted tree admitting the following self-similar 
presentation [8]: 

{01){id,id,a) b = (02) {id, b, id) c ^ (12) {c, id, id), 

where (01), (02) and (12) are transpositions in Sym{3). The associated Schreier 
graphs are self-similar in the sense of |13j . that is, each Sn+i contains three copies 
of Tin glued together by three edges. These graphs can be recursively constructed 
via the following substitutional rules [8]: 
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Remark 3.1. Observe that, for each n > 1, the graph I]„ has three loops, at the 
vertices 0", 1" and 2", labeled by c, b and a, respectively. Moreover, these are the 
only loops in E„. The Ising model will be studied on I]„ considered without loops. 

Let us now proceed to the computation of closed polygons in E,, . Denote by Pn the 
set of closed polygons in E„, and by Ln the set of all subgraphs of I]„ consisting 
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of self-avoiding paths joining the left-most vertex to the right-most vertex in S- 
together with closed polygons having no common edge with the path. 




Two elements of P2 . 




Two elements of L2 ■ 



Each closed polygon in I]„ can be obtained in the following way: either it is a union 
of closed polygons living in the three copies or it contains the three special 

edges joining the three subgraphs isomorphic to The subgraphs of the first 

type can be identified with the elements of the set P^-i, whereas the other ones 
are obtained by joining three elements in i„-i, each one belonging to one of the 
three copies of S„-i, so that they can be identified with elements of the set L^^^. 
This gives 

(2) Pn^PLl]\Ll-l- 

On the other hand, each element in L„ can be described in the following way: if 
it contains a path that does not reach the up-most triangle isomorphic to 
it can be regarded as an element in L/^^^i x P„_i; if it contains a path which goes 
through all three copies of then it is in L^i^i- This gives 

(3) L„= xP„_i)[]lLi, 
from which we deduce 

Proposition 3.2. For each n > 1 the number \Pn \ of closed polygons in the Schreier 
graph S„ of H^^^ is 2^~2-^. 

We are now ready to compute the generating series for closed polygons and the 
partition function of the Ising model on Schreier graphs of H^^\ Denote by T'^jiz) 
the generating function of the set of subgraphs in P„ and by Tn{z) the generating 
function of the set of subgraphs in L„. The equation ([2]) gives 

(4) T';Uz)^{r^_,iz)f + z'Tl_,{z). 

The factor z^ in (jlj is explained by the fact that each term in T^(z) corresponds 
to a set of edges that becomes a closed polygon after adding the three special 
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edges connecting the three copies of I]„-i. We have consequently that the second 
summand is the generating function for the closed polygons containing the three 
special edges. Analogously, from ^ we have 

(5) T„(z) = zT2_,(z)r^'_,(z) + z'Tl_,iz). 

Theorem 3.3. For each n > 1, the partition function of the Ising model on the 
Schreier graph I]„ of the group 

ij(3) 

IS 

Z„ = 2^" • cosh(^J)^^^^^ • r^'(tanh(^J)), 

with 

n 

r:;'w = n ^r' w • (v-n+iw - 1), 
fc=i 

where tpi{z) = and ipk{z) = ip1_i{z) — 3ipk-iiz) + 4, for each k > 2. 

Proof RecaU that Z„ = 2l^«''*(^")l cosh(/3J)l^''9<=''(^")l ■r^'(tanh(/3J)), where rf/(z) 
is the generating function of closed polygons in I]„. In our case we have \Edges(Tin)\ 
= ^^^^ and |yert(I]„)| = 3". 

We know that the generating functions r'^{z) and T„(z) satisfy equations ([4]) 
and ([Sj, and the initial conditions can be easily computed as: 

rf(z) = l + z^ Ti{z) = z^ + z. 

We now show by induction on n that the solutions of the system of equations (|4]) 
and (O are 

frf/(^) = nLi ^r^i^) ■ ii^n+iiz) - 1) 

\t„(z) = .3"-ij^^^^^3"-(,). 

Forn 1, we get (z) = z^ipi{z){ijj2iz)-l) = z'-^ + l andTi(z) = z^^pi^z) = z^+z 
and so the claim is true. Now suppose that the assertion is true for n and let us 
show that it is true for ?i + 1. One gets: 



= iz'"l[^r\z)-{^P„+,{z)-l)\ +z^z'"-'l[^pr\z) 

71 

fc=l 

n+1 



fc=l 

n+1 

fe = l 
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and 



^ / n 



k=l / \ k=l 



n 

fe=l 

n+1 

^3"--l-Q^3"--(,)_ 



k=l 

□ 

Theorem 3.4. The thermodynamic limit lini„_>.oo \veH:f£\\ 

Proof. Since \Edges(T.n)\ = ^ ^2 ^^^^ |^eri(S„)| = 3", the limit reduces to 

log(2) + \ log(cosh(/3J)) + lim 

where tanh(/3J) takes values between and 1. It is straightforward to show, by 
induction, that ipk{z) = '^^^-l 1 for every A: > 1, where fkiz) is a polynomial of 

degree 2*^"^ in z. Hence, the limit lim„^oo '"^^3" becomes 

log (riLl 'pT'\z) ■ (V9„+1 (z) - 

lim 

log(^fe(z)) l0g(<^„+i(z) - 



hm y^^^HV^^ j.^ 



n— J-oo ' 3 71— >oo 3^^ 

Let us show that the series 'J2'kLi '"^^^^f*-^-*-* converges absolutely, and that 
liniji^oo — I — 0. It is not difficult to show by induction that 

2z^'"' < (^fe(z) < 22'-i 

for each fc > 2 and z G [0, 1], so that 

I log(^fc(z))| < max{2'=-i| log(z)|, (2^= - 1)| log(2)|}. 

Note that 1 < (pi{z) < 2 for each < z < 1. Moreover, one can directly verify that 
ip2{z) has a minimum at C2 = 1/4 and (p2{z) < for each z G (0,C2). Let us call 
Cfe the point where ipk{z) has the first minimum. One can prove by induction that 
f'k+iiz) < for each z £ (0,Cfc] and so Ck < Ck+i for every k > 2. In particular, 
(pk{z) satisfies 

-2'=log(2) < log(^fc(z)) < (2^- - l)log(2), 
that gives | log((y3fc(z))| < 2^ log(2) for each k > 2. So we can conclude that 



^ |log(^fc(z))| ^log(2) ^ 2Mog(2) 
3/= - 3 3/= 

k=l k=2 

|log(y„+i(2:)-^"")| ^ 2" log(2) 



< 00. 



Moreover lim„_,oo 1 '°b^'^"+u^^-^ )i < lim„^^ =0. □ 
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3.2. The Sierpinski gasket. In this section we use the high temperature expan- 
sion and counting of closed polygons in order to compute the partition function 
for the Ising model on a sequence of graphs {^ln}n>i converging to the Sierpinski 
gasket. The graphs f2„ are close relatives the Schreier graphs I]„ of the group H^^^ 
considered above. More precisely, one can obtain iln from S„ by contracting the 
edges between copies of S^-i in The graphs ri„ are also self-similar in the 
sense of [13] , as can be seen in the picture. 




Similarly to the case of H^^^ above, define sets P„ and L„. The same recursive 
rules hold, and the total number of closed polygons is again 2^~^ , since the initial 
conditions are the same. 

Let r^'(z) denote the generating function of the subgraphs in P„ and let T„(z) 
denote the generating function of the subgraphs in L„. From relations ([2]) and ^ 
we deduce the following formulae: 

(6) r';l{z) = {r^^_,iz))' + Tl_,{z). 

and 

(7) T„(z) = T2_i(z)rli(z) + Tl_,{z). 

Note that in (jG]) and ([7]) there are no factors z,z'^,z^ occurring in (jlj and ([5]), 
because the special edges connecting elementary triangles have been contracted in 

■ 

Theorem 3.5. For each n > 1, the partition function of the Ising model on the 
n-th Sierpinski graph fi„ is 

Zn = 2^ ■ cosh(^J)3" . r^'(tanh(/3J)), 

with 

n 

where ipi{z) = 1^2(2) = and tpk{z) = i!l_j^{z) ~ 3ipk-iiz) + 4, for each 

k>3. 

Proof Again we shall use the expression Z„ = 2\^ertin„)\ ^Q^Yi{f^jyEdges{n„)\ . 
r^'(tanh(/3J)), where r^'(z) is the generating function of closed polygons in f2„. 
In our case wc have \Edges{nn)\ = 3" and \Vert{iln)\ = ^ 2^ ■ 

We know that the generating functions r^'(z) and Tn{z) satisfy the equations 
(|6]) and ([7]), with the initial conditions 

rf{z)^l + z^ Ti{z) = z^ + z. 
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Let us show by induction on n that the solutions of the system of equations ([6]) and 
dZ]) are 

(r^^iz) = z"^ nLi ^r\z) ■ i^n+iiz) - 1) 
\T„(z) = z^nLiV'f"w- 

For n = 1, we get (2) = z^ipi{z){ip2{z) — I) = z^ + 1 and Ti(z) = zii{ji{z) = 
z'^ + z and so the claim is true. Now suppose that the assertion is true for n and 
let us show that it is true for n + 1. One gets: 

n 

k=l 

n+1 

and 



2 



fc=l / \ fc=l 



k=l 
n 

z"^ Y[^f^\z) ■ ^„+i{z) 



k=l 
n+1 

2 2 Yin {z) 

k=l 



□ 



Remark 3.6. The existence of the thermodynamic limit can be shown in exactly 
the same way as for the Schreier graphs of the Hanoi Towers group. 

3.3. Renormalization approach. Expressions for the partition function of the 
Ising model on the Sierpihski gasket are well known to physicists. A renormalization 
equation for it can be found for example in [5] (see also references therein), and a 
more detailed analysis is given in [2] . Using renormalization, Burioni et al [2] give 
the following recursion for the partition function of the Ising model on the graphs 
n > 1: 

(8) z^+,{y) = zM{y)My)?"~\ 

where y = exp(/3J); f{y) is a substitution defined by 

„8 „,4 , . \ 1/4 



and 

,,4 



{y) = ^-4^[{y' + 3ny'-y' + ^)V/'; 



y3 
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with 

Zi(j/) = 2yV6y-i . 

A similar computation can be performed in the case of the Schreier graph I]„ of the 
group H^'^\ where the self-similarity of the graph allows to compare the partition 
function Zi of the first level with the partition function of level 2, where the sum 
is taken only over the internal spins (T2 1 f 3 j f 5 7 o'e > "'s j "'g (see figure below). The 
resulting recurrence is the same as ([5]), but with 



y f{y) 

and 



2(7TT) 




Remark 3.7. The above recursions for the partition functions for il„'s and S,i's 
can be deduced from our Theorems 13.31 and 13.51 by rewriting the formulae in the 
variable y = exp(/3J) and substituting z ~ tanh(/3J) — (j/^ — + !)■ 



4. Statistics on weighted closed polygons 

This Section is devoted to the study of weighted generating functions of closed 
polygons, i.e., we allow the edges of the graph to have different weights tanh(/3Jij ), 
as in RHS of ([T]). We also take into account the fact that the graphs we consider are 
Schreier graphs of some self-similar group G with respect to a certain generating set 
S', and their edges are therefore labeled by these generators. It is thus natural to 
allow the situations where the energy between two neighbouring spins takes a finite 
number of possible values encoded by the generators S . Logarithmic derivatives of 
the weighted generating function with respect to s G S" give us the mean density of 
s-edges in a random configuration. We can further find the variance and show that 
the limiting distribution is normal. 

4.1. The Schreier graphs of the Grigorchuk's group. Recall from 12.1] that 
the simple edges in I]„ arc always labeled by a. Moreover, the 2-cycles can be 
labeled by the couples of labels (6, c), (6, d) and (c, d). We want to compute the 
weighted generating function of closed polygons, with respect to the weights given 
by the labels a, 6, c, d. Let us set, for each n > 1: 

Xn — I {2-cycles with labels 6, c}| l^n = | {2-cycles with labels 6, c?}| 
Wn = I {2-cycles with labels c, c?}| 
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One can easily check by using self-similar formulae for the generators, that the 
following equations hold: 

In particular, one gets 



with initial conditions Xi 




2. One gets the following values: 



) 2"-2 



- 2 

7 

2"-4 



if 71 = 0(3) 
if 71 = 1(3) 
if n EE 2(3) 



0(3) 
1(3) 
2(3) 



and, consequently. 



Theorem 4.1. For each n > 1, the weighted generating function of closed polygons 
in I]„ is 

'r^'(a,6,c,d) = (1 + bc)^^^^ {I + bd)^ {1 + cdf-^^^ ifn = 0{3) 
r^' (a, 6, c, d) = (1 + be) ^^^^ {I + bd)^ {I + cd) ^^^^ if n = 1(3) 
J^'(a,6,c,d) = (1 + bc)^^^^ {I + bd)^ {I + cd)^^^^ ifn = 2{3) . 

Proposition 4.2. Let w„ be the number of edges labeled w in a random closed 
polygon in E„, where w = a,b,c,d. Denote by ^n,w and tr^ ^ the mean and the 
variance of Wn ■ Then, 

• for each n > 1, On = 0; 

• The means and the variances of the random variables 6„,c„,d„ are given 
in the following table: 





n = 0(3) 


n = 1(3) 


n = 2(3) 




-1(2" -1) 


f(2"-i-l) 


3-2" -5 
14 




i(2" - 1) 


A(2"-i-l) 


3-2" -5 
28 


f-^n,c 


^ 


^(2"-i-l) 


5-2""^-3 





5-2""^-3 
14 


^(2-1-1) 


5-2""^-3 
28 




3.2"-^_.5 



A(2"-i-l) 


f(2«-2 _ 1) 




28 


i(2"-i-l) 





• the random variables 6„,Cn,(i„ are asymptotically normal, as n ^ oo. 

Proof. It is clear that an edge labeled by a never belongs to a closed polygon of 
so that a„ = 0. Let us choose, for instance, n = 0(3). Putting 



V^ib) :-r^'(l,6,l,l) = 2^ 



^(l + b)' 
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we can obtain the mean fin,b and the variance for 6„ by studying the derivatives 
of the function log(r^'(6)). One gets 



^J■n,b 



-(2" 



1) 



— (2" 
28^ 



1) 



Similar computations can be done for c and d. 

In aU cases above we can find exphcitly the moment generating function of the 
corresponding normahzed random variables: a direct computation of its limit for 
n — oo shows that the asymptotic distribution is normal. □ 

4.2. The Schreier graphs of the Basilica group. We also compute the weighted 
generating function of closed polygons for the Basilica group, with respect to the 
weights given by the labels a and b on the edges of its Schreier graph I]„. We use 
here the computations from Proposition [2?2l 

Theorem 4.3. The weighted generating function of closed polygons in the Schreier 
graph I]„ of the Basilica group is 



^nia,b) 



for n > 5 odd and 



n 

fc=i 



1 



n 



fc=i 

2 



l + b^ 



k=l 



for n > A even. 



Proposition 4.4. The means and the variances of the densities 
given in the following table: 





n > 5 odd 


n > 4 even 


t^n,a 


2n-2 




'1 

^n.a 


(n + 1)2"-'^ 


{n + 2)2"-"^ 


tJ-n,b 


2«-i 


2«-i 


^ib 


(n + 3)2"-^ 


{n - 2)2"-^ 



4.3. The Schreier graphs of i/^'^^ Let us denote by Tjj'(a, 6, c) the weighted 
generating function of the subgraphs that belong to the set P„ , defined in Subsec- 
tion [3lT] (the exponent Ir stands for left- right, as self-avoiding paths in P„ join the 
left-most and the right-most vertices of Analogously, we define Tjj'(a, 6, c) and 
T™(a, 6,c), where the exponents lu and ru stand for left-up and right-up, respec- 
tively. By using the self-similar expressions for the generators given in Subsection 
3.1[ we find that these functions satisfy the following system of equations (we omit 
the arguments a, b, c): 



(9) 



n+l - 
^ n+l 
^ n+l 
^ n+l 



^ n ) ' ^ n ^ n n 

T5rT™r^'4-6c(T^")' 



a 

^ ^ n ^ n ^ , . 

\J /y^ X rv-fc ^ /y^ 



ac \ 
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with the initial conditions (a, 6, c) = 1 + abc, T'{{a, 6, c) = ac + b, T']^"(a, 6, c) 
a + be, T™(a, 6, c) = c + afe. 

Proposition 4.5. T/ie mean and the varianee for Wn, with w — a,b,c, are: 

_ 3" - 1 2 _ 3" - 1 

The random variables Wn with w ^ a,b, c are asymptotically normal. 
Proof. If we put a ^ b = 1 , the system ([9]) reduces to 



(10) 



Ti+1 ~ Ti n ri V n J 



with the initial conditions Tf (l,l,c) = T'i''(l,l,c) = T'i"(l,l,c) = T'i"(l,l,c) = 
1 + c. 

One can prove, by induction on n, that the solutions of the system ()10p are 

r:;'(l, 1, c) = T^'^(l, 1, c) = T'„"(l, 1, c) - T™(1, 1, c) - (1 + c)^ for each n. 

By studying the derivatives of the function log(r5^'(l, l,c)) with respect to c, one 
gets: 

3" - 1 2 _ 3" - 1 



Mn,c = cr. 



Symmetry of the labeling of the graph ensures that the same values arise for the 
random variables a„ , 6n . □ 



4.4. The Sierpiiiski graphs. The Sierpihski graphs fi„ being not regular, they 
cannot be realized as Schreier graphs of any group. There exist however a number 
of natural, geometric labelings of edges of r2„ by letters a, 6, c (see [3j). Here we will 
be interested in one particular labeling that is obtained by considering the labeled 
Schreier graph I]„ of the Hanoi Towers group and then by contracting the edges 
connecting copies of I]„_i in E„; and so we call this the "Schreier" labeling of f2„. 

Remark 4.6. The "Schreier" labeling on r2„ can be constructed recursively, as 
follows. Start with the graph VLi in the picture below; then, for each n > 2, the 
graph Vln is defined as the union of three copies of ri„_i. For each one of the 
out-most (corner) vertices of f2„, the corresponding copy of ^n-i is reflected with 
respect to the bisector of the corresponding angle. 
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b c a 5 

Let Tj[(a, &, c), Tj^(a, b, c) and T™(a, &, c) be defined as for the Schreier graphs E„ 
of the Hanoi Towers group in the previous subsection. Then one can easily check 
that these functions satisfy the following system of equations: 



^ n+l - 

^ n+l 

^ n+l 
ylr 



n n n 



with the initial conditions Tf{a, b, c) 
a + be, T™(a, b,c) = c + ab. 

Proceeding as in Subsection 14.31 we find 



n n 

n n n 

^lu'Y'ruT^cl 

n n n y ,^ , 

1 + abc, T5^'(a, 6, c) 



r^/(l, l,c) = TJ:(1, 1, c) = Tjr(l, T™(1, l,c) 

which implies the following 



Proposition 4.7. T/ie mean and the variance for the random variable w„ 
w = a,b,e, for Q„ with the "Schreier" labeling are: 

3«-l ^ gn-l 

l^n.w 



with 



The random variables an,bn, 



2 n.w ^ 

c„ are asymptotically normal. 



It is interesting to compare these computations with those for a different labeling 
of f2„, that we call the "rotation- invariant" labeling of Sierpihski graphs, defined 
recursively as follows. (Compare the construction to the recursive description of 
the "Schreier labeling" in Remark 14.61 ) 

Let be the weighted graph in the following picture. 




a 5 

Then define, for each n > 3, r2„ as the union of three copies of fln-i, rotated by 
kir/S with k = 0, 1,2. 

For rt = 3, one gets the following graph. 
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b n 5- 

It turns out that the weighted generating function of closed polygons is easier to 
compute for 0„ with the " rotation- invariant" labeling, than with the "Schreier" 
labeling. More precisely, we have the following 

Theorem 4.8. For each n > 1, the weighted generating function of closed polygons 
for the graph fin with the "rotation-invariant" labeling is 

2 ^ 

r^'(a, 6, c) = ((a + 6c)(6 + ac))^^ Vf b, c) J| i'f'\a, 6, c) • (V'n+i(a, 6, c) - 1) 



fc=2 



where "01 (^7 c 



l+c 



l+ab 



T, i/'2(a, c) = ■ r 

i(a+bc){b+ac))i {{a+bc){b+ac))2 
aH^c^-aH^c+aH^+4abc+c^^c+l+a^c+b^c r ^ > ^ 

{a+bc){b+ac) ' •' — ' 

'4>k{a,b,c) = i'l_i{a,b,c) - 3ipk^i{a,b,c) +4. 



Proof. Consider the graph fin- For each n > 2, define the sets P„ and L„ as 
in Subsection 13.11 and let rfj(a, 6, c) and T„(a, 6, c) be the associated weighted 
generating functions. By using the symmetry of the labeling, one can check that 
these functions satisfy the following equations 



(11) 



1 r^'(a, 6, c) = (rl,(a, b, c)y + Tl_,{a, b, c) 
\ T„(a, 6, c) - T3_i(a, 6, c) + T2_,(a, 6, c)r$^Li(a, 6, c). 

with the initial conditions 



Aabc — a&c — ab - 



1) 



I r^' = (1 + c)(l + a6)(a262c2 „ ^2^2^ ^2^2 

[T2 = (l + c)(l + a&)(a + 6c)(^ + ac). 

As in the proof of Theorem l3.51 one shows by induction on n that the solutions of 
the system are 

fr^' = ((a + bc){b + ac))^ Vr"(a, c) ^T^a, b, c)-(V„+i(a, b, c) - 1) 

1t„ = ((a + 5c) (5 + ac))^ («, ^ c) UL, i^f-" (a, &, c) 

□ 

Remark 4.9. Although the labels a and are not symmetric to the label c in 
the " rotation- invariant" labeling, computations show that the functions r^'(a, 1,1), 
r^'(l, 6, 1) and r^'(l, 1, c) are the same in this case as in the case of the "Schreier" 
labeling. It follows that the values of the mean and the variance of the random 
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variables a„,6„,Cn remain the same as in the "Schreier" labeling, see Proposition 



4.5. Correspondences via Fisher's Theorem. In |S] M. Fisher proposed a 
method of computation for the partition function of the Ising model on a (finite) 
planar lattice Y by relating it to the partition function of the dimers model (with 
certain weights) on another planar lattice constructed from Y. (The latter par- 
tition function can then be found by computing the corresponding Pfafhan given 
by Kasteleyn's theorem.) This method uses the expression ([T]) for the partition 
function in terms of the generating function of closed polygons in Y. The new 
lattice Y^ is constructed in such a way that Ising polygon configurations on Y are 
in one-to-one correspondence with dimer configurations on Y^. In order to have 
equality of generating functions however, the edges of Y^ should be weighted in 
such a way that the edges coming from Y have the same weight tanh(/3Ji j ) as in 
the RHS of (ID), and other edges have weight 1. 

Applying Fisher's construction to Sierpiiiski graphs, one concludes easily that if 
Y = fin for some n > 1, then Y"^ = S„+i, the (n-l-l)-st Schreier graph of the Hanoi 
Towers group i?'-^-' with three corner vertices deleted. Note that the corner vertices 
are the only vertices in E„ with loops attached to them, and so it is anyway natural 
to forget about them when counting dimer coverings. The construction consists in 
applying to Y the following substitutions, where edges labeled by e in are in 
bijection with edges in Y, and should be assigned weight tanh{/3Jij). Other edges 
should be assigned weight 1. 




The correspondence between closed polygons in Y and dimer coverings of Y^ is as 
follows: if an edge in Y belongs to a closed polygon, then the corresponding e-edge 
in does not belong to the dimer covering of Y^ associated with that closed 
polygon, and vice versa. 

The following pictures give an example of a closed polygon in il2 and of the 
associated dimer covering of E3 : 
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If, for a certain n > 1, the Sierpinski graph fi„ is considered with the "Schreier" 
labeling, then the labeling of the graph S„+i given by Fisher's construction will be 
a restriction of the usual Schreier labeling of More precisely, only the edges 

that connect copies of but not copies of S,i will be labeled (other edges have 

weight 1), and the labels are the same as in the standard labeling of E„+i as a 
Schreier graph of the group TJ^'^^. The following picture represents S3 as with 
F = O2 with the "Schreier" labeling. 

i/\i i/\i 




Remark 4.10. One can wonder what Fisher's construction gives for {S„}„>i. It 
turns out that if F = S„, the rt-th Schreier graph of H'^^\ then = the 
same as for Y = f2„, the n-th Sierpinski graph . 
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